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Abstract. We present a detailed phenomenological study of the prospects of open charm physics at the 
future pp and pp facilities GSI-FAIR and J-PARC, respectively. In particular, we concentrate on differential 
cross sections and the charge and longitudinal double-spin asymmetries at next-to-leading order accuracy. 
Theoretical uncertainties for the proposed observables are estimated by varying the charm quark mass and 
the renormalization and factorization scales. 
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In recent years, the study of heavy flavors at colliders has 
become a versatile tool to probe different aspects of Quan- 
tum Chromodynamics (QCD), ranging from heavy flavor 
parton densities and the hadronization of heavy quarks 
into heavy mesons or baryons to the dynamics of QCD 
hard-scattering [T]. Technological advances, like the intro- 
duction of micro-vertex detectors, allow for much better 
tagging of the produced heavy quarks and hence more 
precise measurements. An equally significant improvement 
of the theoretical foundations for heavy flavor production 
has been achieved, in particular, novel proposals to over- 
come limitations of fixed-order calculations and to model 
the hadronization of heavy quarks more reliably PQ . As a 
consequence, a discrepancy between data and theory for 
bottom production, suggested by various experiments, has 
been reduced to the point that it no longer appears sig- 
nificant [2]. 

A solid understanding of the theoretical framework 
for heavy quark production and its phenomenological in- 
puts is of utmost importance for several key measurements 
both ongoing or taking place in the very near future [1] . At 
the LHC, heavy flavor production by genuine QCD pro- 
cesses is an important background to searches for Higgs 
bosons within the Standard Model and beyond, e.g., to the 
decay H — > bb, as well as for new physics. At the Rela- 
tivists Heavy Ion Collider (RHIC) one wants to establish 
the existence and to investigate the properties of a new 
state of matter, the quark-gluon plasma (QGP). Several 
signatures related to heavy flavors have been proposed, in 
particular, how their production and decays are modified 
by the presence of a QGP pj]. Here, heavy flavor produc- 
tion in proton-proton collisions will act as a benchmark for 
possible modification in nucleus-nucleus collisions at the 
same energy [3,4 . In addition, RHIC is also capable to 
collide longitudinally polarized protons at center-of-mass 



system (c.m.s.) energies of up to a/5 = 500 GeV. Besides 
its intrinsic interest as an important test of the dynamics 
of spin-dependent hard-scattering in QCD, these measure- 
ments are likely to further our understanding of how glu- 
ons contribute to the proton spin, a major goal of hadron 
physics [5]. 

In terms of perturbative QCD (pQCD), relevant fixed- 
order calculations have been in place for quite some time 
now, and next-to- leading order (NLO) accuracy is the 
state-of-the-art throughout. NLO results, which keep the 
full dependence on the heavy quark mass ttiq , for unpolar- 
ized hadron-hadron ("hadroproduction"), photon-hadron 
("photoproduction"), and electron-hadron ( "electroproduc- 
tion") collisions can be found in [6], [7], and [8], respec- 
tively. Calculations for longitudinally polarized hadro- [9] 
and photoproduction [lOj are more recent achievements. 
In each case one exploits the fact that uiq acts as an effec- 
tive infrared cut-off for collinear singularities, which allows 
to compute total heavy quark yields as a perturbative se- 
ries in the strong coupling a s evaluated at a hard scale 
of 0(jtiq). For inclusive transverse-momentum pt distri- 
butions, complications arise once pt/tuq 3> 1, and large 
logarithms in this ratio have to be resummed to all orders 
in a s to improve the convergence of the perturbative series 

mi. 

In this paper, we will focus on a detailed phenomeno- 
logical study of the prospects of open charm production 
in pp and pp collisions at the future GSI-FAIR l 12] and 
J-PARC P3] facilities, respectively. For GSI-FAIR we con- 
sider the so-called "collider option" as proposed by the 
PAX collaboration [14], using a 15 GeV anti-proton and 
a 3.5 GeV proton beam, with the plan of having both 
beams also longitudinally or transversely polarized. The 
main goal of the PAX experiment would be a determina- 
tion of the so far unknown "transversity" parton densi- 
ties through the transversely polarized Drell-Yan process. 
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This measurement would gain substantially from a polar- 
ized anti-proton beam due to the dominance of the low- 
est order (LO) quark- antiquark annihilation channel. The 
J-PARC facility in Japan is currently being completed. 
Here we consider the collision of a 50 GeV proton beam 
with a fixed, solid-state target. The possibility of having 
both beam and target polarized is a conceivable option 
for future upgrades currently being scrutinized. At GSI- 
FAIR (J-PARC) the available pp (pp) c.m.s. energy will 
be y/S ~ 14.5 (10) GeV. For both experiments the details 
of the detector and the acceptance are not yet finalized. 
We will make some reasonable assumptions as stated be- 
low in Sec. 13.11 We will demonstrate that measurements 
of open charm at these facilities have the potential to fur- 
ther our understanding of the underlying QCD dynamics 
at moderate c.m.s. energies so far little explored. 

First, one has to determine, of course, to what extent 
perturbative QCD is applicable at c.m.s. energies of about 
10 to 15 GeV. NLO calculations for single-inclusive pr- 
spectra of pions or photons are known to seriously under- 
shoot data even at somewhat higher c.m.s. energies [To] . 
The mass of the charm quark already sets a hard scale 
of 0(1 GeV), which may facilitate the use of perturbative 
QCD even for small or vanishing transverse momenta pr 
of the observed charm quark. Also, since px ~ C( to q)i 
we do not have to worry about potentially large loga- 
rithms lnp<r/mQ present at large c.m.s. energies |llj . how- 
ever, partonic threshold effects may become important 
and perhaps need to be resummed to all orders. In case 
of px-differential pion spectra, it was shown that thresh- 
old resummations can lead to a much improved agree- 
ment between theory and experiment at low c.m.s. ener- 
gies [15] . To study the possible relevance of resummations 
at GSI-FAIR and J-PARC, we provide total and differ- 
ential charm yields for both pp and pp collisions at NLO 
accuracy of QCD, including detailed discussions of theo- 
retical uncertainties due to variations of the charm quark 
mass or the renormalization and factorization scales. We 
note that the computation of threshold resummations for 
charm production at GSI-FAIR and J-PARC is far beyond 
the scope of this paper, which aims at a first exploratory 
study of charm physics at comparatively low c.m.s. ener- 
gies. 

Other interesting, though experimentally more chal- 
lenging observables related to open charm production can 
be considered at GSI-FAIR and J-PARC. One is the so 
called "charge asymmetry" which describes the difference 
of cross sections for producing a heavy quark Q or a heavy 
anti-quark Q at a certain point in phase-space: 



da® - do® 

da® + daQ 



(1) 



This asymmetry probes a subset of NLO radiative correc- 
tions and vanishes at the LO approximation. This feature 
makes it an important test of QCD hard scattering dy- 
namics. The Abelian (QED) part of Ac is also known as 
forward-backward asymmetry |16| and is caused by the 
interference of states with different C-parity. It should be 
mentioned that routinely used event generators based on 



LO matrix elements cannot predict this interesting effect. 
In the unpolarized case the charge asymmetry was first 
mentioned in [T7l[6] and later studied qualitatively in [TBI 
[19], mainly for top production at the TeVatron and the 
LHC. A first measurement of Ac for top production has 
been recently reported by the CDF collaboration [20j . We 
estimate the size of this effect for GSI-FAIR and J-PARC 
and, for the first time, compute the corresponding charge 
asymmetry also for polarized hadroproduction. We will 
give more details in Sec. 12.21 

Provided that longitudinally polarized beams and tar- 
gets will be available at GSI-FAIR and/or J-PARC, stud- 
ies of the double-spin asymmetry 



All = 



dAa 
da 



(2) 



could provide unique insight in the distributions describ- 
ing the polarization of quarks and gluons in the nucleon 
at medium-to-large momentum fractions x. In @, dAa 
denotes the spin-dependent cross section, which can be 
obtained by taking the difference of measurements with 
the (anti-)proton spins aligned and anti-aligned. The un- 
polarized cross section da is determined by the sum of 
both measurements. We will explore the sensitivity of All 
to different assumptions about polarized parton densities, 
in particular, the elusive gluon polarization. As for the 
unpolarized cross sections, we discuss theoretical uncer- 
tainties related to variations of the charm quark mass and 
the renormalization and factorization scales. 

The paper is organized as follows: in Sec. 2 we briefly 
review the technical framework for hadroproduction of 
heavy flavors, set our notation, and introduce in some 
detail the charge asymmetry. Section 3 is devoted to nu- 
merical studies. First we define the phcnomenological in- 
puts and experimental cuts we assume for GSI-FAIR and 
J-PARC. Next, we give results for total and differential 
unpolarized charm cross sections and discuss their theo- 
retical uncertainties. These results serve as a future bench- 
mark to study the applicability of pQCD for charm pro- 
duction at yS ~ 10 GeV. Then we present expectations 
for the polarized cross sections and charge and double spin 
asymmetries. We briefly summarize the main results in 
Sec. 4. 



2 Technical Framework 

The pQCD framework for single-inclusive heavy flavor 
production at NLO accuracy in both unpolarized and lon- 
gitudinally polarized hadron-hadron collisions has been 
discussed in detail in Refs. [6] and [9j, respectively. We 
can restrict ourselves to a brief summary of the aspects 
with particular relevance to our analysis. 



2.1 The single-inclusive cross section 

We are interested in the hadroproduction of a heavy quark 
Q [anti-quark Q] with mass vtlq and four-momentum Pq 
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Hi(Pi) J ff 2 (P 2 ) - Q(Pq)[Q(Pq)] + X. (3) 

X includes all other final-state particles such that (J3j> is 
inclusive with respect to the detected heavy quark. Pi 2 
denote the momenta of the colliding hadrons H12 and 
S = (Pi + P2) 2 the c.m.s. energy squared. 

More specifically, Pq can be parameterized in terms 
of the transverse momentum px, the rapidity y, and the 
azimuthal angle <fi of the observed heavy quark Q, 



Pq = (rriT cosh y, px sin <f>, px cos <p, mx sinh 



y, 



(4) 



where the transverse mass is defined as mx = yPx + to q 

and y = \ ln[(P + P Z )/(E - P z )]. Applying the factoriza- 
tion theorem, the unpolarized differential cross section for 
([3]) can be written schematically as 



drrvLdy 



ab 



d 2 a ab ^Qx (s 1 , x 2 , Pi , P-i , Pq , \if , \i r ) 
dm T dy 



(5) 



with the symbol ® indicating a convolution. X\ and Xi are 
the fractions of Pi and P2 taken by the partons f^ 1 and 
f b H2 with flavor a and b, respectively. The sum in Eq. |5]) 
is to be performed over all contributing partonic channels 
ab — > QX with da a b-,Qx the associated partonic cross 
sections. 

The factorization of the cross section (J5J) into non- 
perturbative parton densities f a .b and short-distance cross 
sections requires the introduction of factorization and re- 
normalization scales /i/ and \i r , respectively. These scales 
are essentially arbitrary and usually chosen to be of the 
order of a hard momentum transfer characterizing the pro- 
cess under consideration. Any residual dependence of the 
right-hand-side (r.h.s.) of Eq. (JSJ) on the actual choice for 
fif^ r represents an important part of the uncertainties in 
the theoretical description of ([3]). 

The basic framework outlined above carries over to 
the case of polarized hadron-hadron collisions as well. To 
obtain the spin-dependent cross section d 2 Aa / dm T dy en- 
tering the experimentally relevant spin asymmetry All 
defined in Eq. ^ , the parton densities and hard scattering 
cross sections on the r.h.s. of Eq. §5§ have to be replaced 
by their polarized counterparts Af a j, and dA<j a b^QX, re- 
spectively. 

Knowledge of higher-order corrections in the perturba- 
tive expansion of the partonic cross sections da a b^QX and 
dAa a b^Qx is generally indispensable. On the one hand, 
in hadronic scattering they are often sizable, and, on the 
other hand, they are expected to reduce the artificial de- 
pendence on the choice of fXf >r - Also, higher-order correc- 
tions affect unpolarized and polarized cross sections dif- 
ferently and hence do not cancel in the ratio All- 

At 0(a 2 ), the LO approximation, the hadroproduc- 
tion of heavy quarks proceeds through only two partonic 
channels 

99 —> QX and qq — > QX (6) 



where X = Q. The 0(a s ) radiative corrections to ([6|) com- 
prise additional real gluon emission, X = Qg, as well as 
one-loop (virtual) contributions. In addition, a new type 
of subprocess, gluon-[anti]quark-scattering, gq[q] — > QX 
has to be considered at NLO. A detailed account of the 
calculation of the relevant matrix elements, the required 
loop and phase-space integrations, and the cancellation of 
singularities is given in Refs. [B] and [5]- 

Since tuq acts as an effective cut-off for collinear sin- 
gularities, also total heavy quark yields are amenable to 
pQCD. They are obtained by integrating ([5]) over the en- 
tire phase-space using 



S/4 



dm T 



cosh- 1 Vs/(2 niT ) 



•y — COSll 

!J ln T^f / .5/(4cosh 2 y) 

dy 



dy 



2 111 1-/3 



dm T , 



(7) 



where j3 = y 1 — 4itlq/S. Alternatively, one can first de- 
rive the so-called "scaling functions", [Z\]er a h(s, m^), as a 
function of the partonic c.m.s. energy s by integrating the 
partonic cross sections d[A]a a b^Qx *>.!) . which then in 
turn have to be combined with the appropriate combina- 
tion of parton densities. We note that the experimental 
determination of the total cross section always involves 
some extrapolation beyond the accessible ranges in px and 
y and is therefore less reliable and useful than differential 
rates in testing pQCD predictions. In our numerical stud- 
ies for J-PARC and GSI-FAIR we will therefore mainly 
focus on differential cross sections, which are also relevant 
for the charge asymmetry. 



2.2 The charge asymmetry 

At the LO approximation the processes © relevant for 
heavy flavor production do not discriminate between a 
produced heavy quark Q and a heavy anti-quark Q. Hence, 
at any given point (px, y) in phase-space the yields ([5]) for 
Q and Q are identical. 

Radiative corrections change this picture and give rise 
to the charge asymmetry Ac , defined in Eq. (TT]) . Any mea- 
surement of Ac will directly probe and perhaps improve 
our understanding of QCD dynamics beyond the LO. So 
far this higher-order effect has received relatively little at- 
tention [17,18,15] but was recently measured for the first 
time in case of top production at CDF [20], the asym- 
metry exceeding the few percent predicted in [18] with 
large experimental uncertainties though. We will explore 
the prospects of accessing Ac in pp and pp collisions at J- 
PARC and GSI-FAIR, respectively. Due to the relatively 
low c.m.s. energies available, we have to limit ourselves 
to charm quark production. Given that longitudinally po- 
larized beams and targets are a viable option, we also 
explore the polarized charge asymmetry, which we define 
as in Eq. (fT]) with the cross sections da®^ replaced by 
their polarized counterparts dAa®^. 
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Fig. 1. Sample qq cut diagram contributing to Ac 
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Fig. 2. Sample gq cut diagram contributing to Ac- 



At NLO the charge asymmetry receives only a very 
limited number of contributions. Instead of making use of 
the results available in the literature [BJIH] > we chose to re- 
calculate the numerator of An from scratch in both the un- 
polarized and polarized caso The origin of Ac resides in 
the Abelian (QED) part of the higher order processes and 
relates to the interference between amplitudes which are 
relatively odd under the exchange of Q and Q. The gluon- 
gluon fusion process in ([6]) is evidently charge symmetric 
also beyond the LO and does not contribute to Ac- Also 
non- Abelian amplitudes involving the triple-gluon vertex 
depend only on Q + Q and lead to contributions which are 
symmetric under exchanging Q and Q. Upon close exam- 
ination of the remaining amplitudes for the qq and gq [gq] 
initiated subprocesses at NLO, one finds that only such 
"cut diagrams" contribute which have three vertices on 
both the heavy and light quark line. In case of the virtual 
loop corrections to the LO qq process, only the interference 
of the box diagrams with the Born amplitude contributes. 
Examples of cut diagrams relevant for the computation of 
Ac and A Ac are depicted in Figs. [1] and [U 

This observation can be readily understood. To this 
end, let us write the partonic subprocess cross section 
d<7 a b—>Qx as the sum of the interference contributions of 
all contributing amplitudes labeled by the indices i,j 



ib^QX 



(8) 



For simplicity, pre-factors such as the spin and color av- 
erages are mapped into K a \y. Phase-space integration for 
all unobserved partons is implicitly understood in (JSJ) . Ex- 
pressions similar to (jHJ) hold for polarized partonic cross 
sections dAa a b^Qx as well as for ab — > QX. For the nu- 
merator of Ac we have to examine the difference 



Aij = MiM* 



ab^QX 



M,M' 



ab^QX 



(9) 



1 The lengthy expressions are available upon request from 
the authors. 



for all i, j. By interchanging Q and Q, the Dirac struc- 
ture relevant for MiM* changes sign for an odd number 
of propagators in the heavy quark trace and otherwise 
remains the same. The color structure is invariant, ex- 
cept for Tr[T fc T ; T m ] = (d Hm + i/jw m )/4 which contains a 
symmetric and an anti-symmetric piece, dkim = d m ik and 
fkhn = -fmik, respectively. 

Combining everything, only those cut diagrams with 
three vertices on both the heavy and the light quark line, 
cf. Figs. [T] and [21 contribute to Ac, and one finds [18] 



A^ = -(dkim) 2 MiM* 



ib^QX 



(10) 



with (dki m ) 2 — 40/3. MiM* denotes the interference of 
the two amplitudes with the QCD color structure taken 
aside, and is the same (up to pre-factors) as for corre- 
sponding QED processes [IS] , e.g, e + e~ — > Let 
us mention that there is a similar effect in the QCD scale 
evolution of parton densities at next-to-next-to-leading or- 
der (NNLO) also proportional to (dkim) 2 and leading to a 
strange quark asymmetry s(x) ^ s(x) [21] . 

All contributions to the numerator of Ac at O(a^) are 
free of ultraviolet as well as collinear singularities as a con- 
sequence of the symmetry of the LO processes ([6J under 
exchanging Q and Q. Infrared (IR) singularities appear 
in both real gluon emission and virtual loop corrections 
to the LO qq process and cancel in the sum. Effectively 
this implies that the NLO matrix elements only contain 
the charge asymmetry at LO approximation. As in [6,9 
we use dimensional regularization to deal with IR poles in 
intermediate steps of the calculation. In [18] a small gluon 
energy E 9 cut was used to cut off IR singularities. While 
the charge asymmetry also appears in the limit rug. — > 0, 
it vanishes for the total heavy quark cross section as a 
consequence of charge conjugation invariance. 

In Sec. 13.31 we give some quantitative predictions for 
A c and AA C at J-PARC and GSI-FAIR. Due to the low 
c.m.s. energies, gluons are much less abundant than at 
high-energy colliders. Since gluon- gluon fusion only con- 
tributes to the denominator of Ac, the prospects for study- 
ing Ac can be in fact more favorable at GSI-FAIR or J- 
PARC than at colliders like the TcVatron or the LHC. 
Here, gluon-gluon fusion is by far the dominant mech- 
anism for heavy quark production, and one has to find 
suitable corners of phase-space to make Ac experimen- 
tally accessible. One important possibility is top quark 
production .18,19,20 , which receives important contribu- 
tions from qq annihilation also at large S thanks to the 
sizable tuq. 



3 Phenomenological Applications 

Before turning to the prospects of heavy flavor physics at 
GSI-FAIR and J-PARC, we quickly review the informa- 
tion on charm production in pp collisions at c.m.s. ener- 
gies below = 50 GeV gathered so far. Figure [3] shows 
the available data [22] compared to our calculations at 
NLO accuracy using different values of the charm quark 
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Fig. 3. Experimental results [22] for the total charm produc- 
tion cross section at fixed-target energies compared to NLO 
pQCD calculations for three different values of the charm quark 
mass m c . In each case, the shaded band indicates the theo- 
retical uncertainties from varying /i = fi r = fif in the range 
m c < /i < 2m c . 

mass and scales fi = [if = [j, r in Eq. ([5]). As can be seen, 
uncertainties from small variations of m c are as impor- 
tant as scale ambiguities, and both combined can lead to 
almost an order of magnitude change in the total charm 
yield at y/S ~ 10 + 20 GeV. From Fig. [3] one can also infer 
that the theoretical uncertainties become somewhat less 
pronounced with increasing c.m.s. energy. 

Although most experimental results can be described 
with the same choice of m c and scale /x, there is a clear 
need for further precise measurements, in particular closer 
to threshold, below the result of NA32, which does not 
line up so well with other experiments. GSFFAIR and J- 
PARC, to which we shall turn now, can explore this energy 
range in the future. 



3.1 Phenomenological inputs and experimental cuts 

The part of the future GSI-FAIR accelerator complex ame- 
nable to pQCD studies is an asymmetric proton-antiproton 
collider option proposed by the PAX collaboration [TJ] 
with maximum beam energies for protons and antiprotons 
of 3.5 GeV and 15 GeV, respectively, resulting in a c.m.s. 
energy of about y/S = 14.5 GeV. Studies of methods to 
polarize both beams either longitudinally or transversely 
are currently pursued [14]. The PAX detector will have 
nearly full azimuthal acceptance and a polar angle cover- 
age from 5 to 130 degrees is envisioned [T4"Il2"3"] . 

The proton accelerator at J-PARC, which will reach up 
to 50 GeV beam energy, is currently under construction, 
and the hadron physics programme will commence in the 
near future at a c.m.s. energy of about \/S = 10 GeV with 



both beam and target being unpolarized. We also consider 
longitudinally polarized collisions, which are a conceivable 
upgrade of the J-PARC facility in the future [13]. Since 
not very much is known about the experimental set-up at 
this stage, we assume a forward spectrometer geometry 
with a 200mrad acceptance, similar to the one used by 
the COMPASS experiment at CERN. 

Since the details of charm detection in experiment are 
not yet available, we will perform all our calculations on 
the charm quark level, i.e., we do not attempt to model 
the hadronization of charm quarks into charmed mesons 
and their subsequent decays. For the PAX experiment, 
however, identification of open charm events most likely 
proceeds through the detection of decay muons for which 
a lower limit on their momentum of p M = 1 GeV is re- 
quired 14.23J. Therefore we impose a similar cut on the 
laboratory momentum of the primary charm quark in all 
our calculations for GSI-FAIR. 

In all unpolarized calculations at LO and NLO accu- 
racy we use the CTEQ6L1 and CTEQ6M parton distri- 
bution functions [24] and the corresponding LO and NLO 
values for the strong coupling, respectively. Unless stated 
otherwise, the GRSV "standard" set of helicity-dependent 
parton densities [25] is used in the computation of polar- 
ized cross sections and the charge and spin asymmetries. 
Since m c = 1.35 GeV provides a good description of most 
of the data shown in Fig. [3] we make this our default 
choice for the charm quark mass. 

3.2 Expectations for charm production cross sections 

Figure [4] shows our expectations for the unpolarized charm 
production cross section ([5]) at GSI-FAIR at LO and NLO 
accuracy, integrated over transverse momentum px > p™ m 
and the angular acceptance of 5° < C < 130° for PAX. 
The shaded bands indicate the theoretical uncertainties 
when the factorization and renormalization scales are var- 
ied simultaneously in the range m™ m < /if — [i r < 2m™ m . 
Also shown in Fig. [4] are the "K-factors" , the ratio of 
the NLO and LO cross sections, for two choices of scales 
fj,f = fj, r and the fractional contributions of the different 
partonic channels afj to the NLO cross section. 

Besides the sizable dependence on the scales there 
is also a similar uncertainty due to the choice of m c in the 
region p™ ln < 1 GeV, in line with the observations for the 
total charm yields in Fig. [3] For p™ m > 2 GeV, however, 
varying m c in the range 1.25 GeV < m c < 1.45 GeV has a 
negligible impact on the cross section shown in the upper 
panel of Fig. [4] It is worth to notice that there is only 
a rather marginal reduction in the scale ambiguity when 
going from the LO to the NLO approximation. This is 
not unexpected for experiments with limited c.m.s. ener- 
gies, and similar observations have been made for single- 
inclusive hadron production |26] . From the lower panel of 
Fig. S] one can infer that the quark- antiquark annihilation 
subprocess is the most important contribution to the cross 
section. This can be readily understood since quarks and 
antiquarks are both "valence" partons in the proton and 
antiproton, respectively, and from the fact that one probes 
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Fig. 4. upper panel: LO and NLO unpolarized charm cross 
section at GSI-FAIR, integrated over pr > p™ ln an d the an- 
gular acceptance 5° < 9 C < 130°, using m c = 1.35 GeV. 
The shaded bands indicate the uncertainties from varying 
H — fi r = fif in the range mf ln < fj, < 2m™ m ; middle panel: ra- 
tio of the NLO and LO cross sections ( "K- factor" ) ; lower panel: 
fractional contribution of the different partonic channels of,- to 
the NLO cross section a c for ^ = \/2m™ ln . 



fairly large momentum fractions x\ 2 ^ 0.1. The genuine 
NLO quark-gluon subprocess is negligible in the entire px 
range shown. 

The corresponding results for J-PARC are summarized 
in Fig. [3] There are striking differences compared to the 
result for GSI-FAIR shown in the previous figure. The de- 
pendence of the charm cross section on unphysical scales 
fif^ r is even larger here and does not improve when NLO 
corrections are included. This can be taken as a strong in- 
dication that higher order terms in the perturbative series 
are very important. Hopefully, a resummation of the lead- 
ing terms to all orders in a s will tame the scale ambigui- 
ties. Secondly, the size of the NLO corrections compared to 
the LO term, displayed in the middle panel of Fig.O seems 
to be beyond control. Most of the pathological behavior 
of the if-factor at large p^ m can be attributed to the dif- 
ferences in the LO and NLO gluon distributions at large 
xi.2, where they are basically unconstrained by data [24.. 
If one uses NLO parton densities in the calculation of the 
LO cross section, the if-factor does not show such a sharp 
rise, though it remains large (dotted curve in Fig. [5]). De- 
spite the rather large x\ t % values probed at J-PARC, and 
contrary to what happens in pp collisions at GSI-FAIR, 
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Fig. 5. Same as in Figure [4] but now for J-PARC kinematics. 



the gluon-gluon fusion subprocess is by far the dominant 
mechanism to produce the charm quark. We also note that 
varying m c in the range 1.25 GeV < m c < 1.45 GeV has 
a somewhat bigger impact on the cross sections displayed 
in the upper panel of Fig.[S]than in case of GSI-FAIR due 
to the smaller c.m.s. energy. 

Figures [3] and O already demonstrate the potential of 
future low energy pp and pp experiments in further con- 
straining quark and gluon distributions, respectively, in 
the medium-to-large x region, difficult to access at high 
energy colliders. They also show, however, that the pQCD 
framework can not be taken for granted in this energy 
regime, and its applicability has to be carefully scruti- 
nized first by comparing the theoretical expectations with 
data. 

Figure [6] shows the corresponding longitudinally polar- 
ized cross sections at NLO accuracy for GSI-FAIR and J- 
PARC, respectively. Instead of giving also the LO results 
[the if-factors are smaller than in the unpolarized case, 
between 1.5 and 2 (2 and 4) for GSI-FAIR (J-PARC)], 
we chose to display the dependence of the polarized cross 
sections on m c . Since J-PARC can cover only a smaller 
range in pt than GSI-FAIR due to the smaller y/S, the 
mass effects are more important here. The fractional con- 
tributions of the different subprocesses are very similar to 
those shown in Fig. [4] for GSI-FAIR and strongly depen- 
dent on the size of the polarized gluon distribution in case 
of J-PARC, as can be expected already from the lower 
panel of Fig. [5l We note that at J-PARC kinematics, the 



J. Riedl, A. Schafer, M. Stratmann: Prospects of Open Charm Production at GSI-FAIR and J-PARC 



7 




i i i i i i i i i ' i i i i ' i i i i ' i i i i 



5 



mm vr -, -i n 

P T [GeV] 



2 3 

min Tr ~, -ttt 

P T [GeV] 



Fig. 6. Expectations for the longitudinally polarized cross section for charm production at GSI-FAIR (left panel) and J-PARC 
(right panel), using the GRSV "standard" set |25| . Note that the results on the left hand side are for —Aa. The scale uncertainty 
is shown for three different choices of m c , varying y, = y r = (if in the range m™ m < (i < 2m™ m . 



polarized cross section exhibits a node at p™ ln ~ 1 GeV 
if the GRSV " standard" parton densities are used in the 
calculation. 



3.3 The unpolarized and polarized charge asymmetry 

We now turn to a detailed discussion of the charge asym- 
metry Aq defined in Eq. (JXJ) and Sec. 12.21 We show expec- 
tations for the size of the effect and discuss the theoretical 
uncertainties due to variations of fif <r and m c . All results 
are presented as a function of the c.m.s. rapidity y of the 
heavy (anti-)quark, which is related to the rapidity in the 
laboratory frame yi a b by a simple additive boost. Positive 
rapidities refer to the direction of the anti-proton and pro- 
ton beam at GSI-FAIR and J-PARC, respectively. Recall 
that rapidity y and pseudo-rapidity rj are not the same for 
massive particles. The relation between y and the scatter- 
ing angle 9 C of the heavy quark depends both on px and 

77V 
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■ tanh y 




1 — tanh y 

while cos0 c = tanh r). Therefore, angular cuts imposed 
upon the heavy quarks by the experiments do not trans- 
late in simple, px independent cuts for rapidity-dependent 
differential cross sections. 

The upper panels of Figs. [7] and [H] show our expecta- 
tions for the unpolarized charge asymmetry Ac at 0{or s ) 
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Fig. 7. The unpolarized charge asymmetry Ac (upper panel) 
and the NLO c.m.s. rapidity-dependent differential charm plus 
anti-charm cross section da c /dy+da c jdy (lower panel) for GSI- 
FAIR. The scale uncertainty is shown for three different choices 
of m c , varying (i = (i r = (if in the range m c < fj, < 2m c . 
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Fig. 8. Same as in Fig. [7] but for J-PARC kinematics. Fig. 9. Same as in the upper panels of Figs. [7] and [8] but now 

for the polarized charge asymmetry AAc- 



for charm quarks at GSI-FAIR and J-PARC, respectively, 
using the phenomenological inputs and experimental ac- 
ceptance cuts specified in Sec. 13.11 The results for Ac are 
largely independent of the choice for the charm quark mass 
m c , in contrast to the sizable mass dependence observed 
for the NLO c.m.s. rapidity-dependent differential cross 
sections for the sum of charm and anti-charm production 
(lower panels), which enters in the denominator of Ac in 
Eq. Q. 

Note that the c.m.s. rapidity y and the experimen- 
tally relevant rapidity yi ab in the laboratory frame are 
simply related by yi ab = y - 0.737 and y lab = y + 2.334 
for GSI-FAIR and J-PARC, respectively. The scale nj >T 
dependence partially cancels out in Ac, as can be seen by 
comparing the upper and lower panels of Figs. [7J and [HI A 
residual dependence on /x $ T is not surprising since the nu- 
merator of Ac at 0(al) is effectively a LO approximation. 
It vanishes at 0(a 2 s ) and is free of collinear singularities. 
Also note that the sharp drop of da c /dy + da c /dy in the 
lower panel of Fig. [7] is due to the cut imposed on the 
momentum of the heavy (anti-)quark, pi a b > 1 GeV. 

Figure [9] shows our results for the longitudinally po- 
larized charge asymmetry AAc defined as in Eq. |T]), but 
with all cross sections da replaced by their helicity depen- 
dent counterparts dAa. Again, the scale uncertainty is 
shown for three different choices of m c , varying fi = fj, r = 
\if in the range m c < fi < 2m c . All results are obtained 
with the GRSV "standard" set [JS] of spin-dependent par- 
ton densities. As in the unpolarized case, the mass depen- 
dence largely drops out in AAc, in particular for GSI- 
FAIR, but a residual scale fXt r ambiguity remains. 

Since the mass and scale dependence of dAa c /dy + 
dAa c I dy is qualitatively very similar to the corresponding 
unpolarized cases shown in the lower panels of Figs. [Jj and 
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Fig. 10. Numerator (upper panel) and denominator (lower 
panel) of the polarized charge asymmetry AAc at J-PARC, 
calculated with different sets of polarized parton densities. 



[8j we refrain from giving these results here. Instead, we 
shall discuss the significant dependence of both numerator 
and denominator of AA C on the choice of a particular set 
of polarized parton densities. 

Figure [TU] shows the numerator and the denominator 
of the polarized charge asymmetry AAc a t J-PARC, ob- 



J. Riedl, A. Schafer, M. Stratmann: Prospects of Open Charm Production at GSI-FAIR and J-PARC 



9 



-0.1 



-0.2 



-0.3 - 



-0.4 - 



-0.5 




A 



LL 



m c = 1.45 GeV -*m 
m c = 1.35 GeV 
m c = 1.25 GeV 
i i i i I i i i i I i i i i I i i 



1 




; 1 J-PARC 



5 



mm vt ~, -i n 

P T [GeV] 



0.06 



- 0.04 



- 0.02 



-0.02 



- -0.04 



mm r/ -, -i J-, 

P T [GeV] 



Fig. 11. The integrated double-spin asymmetry All as function of p™ 111 at NLO accuracy for GSI-FAIR (left panel) and J- 
PARC (right panel) using the GRSV "standard" distributions and the same experimental cuts as before. The scale uncertainty 
is shown for three different choices of m c , varying fi = /i r = fif in the range m™ m < /J, < 2m™ m . 



tained with different sets of polarized parton densities. 
Apart from our default set, GRSV "standard" , we also use 
the sets of AAC [27] and DNS [55]. The latter is based on 
an analysis using also data from polarized semi-inclusive 
deep-inelastic scattering. Two different choices of parton- 
to-hadron fragmentation functions have been made in the 
DNS analysis, and the two resulting sets, labeled as DNS 
(KRE) and DNS (KKP) in Fig. [TO] differ mainly in the 
sea quark content, in particular, Au, which has oppo- 
site sign in both sets. The positive polarization of Au in 
DNS (KKP), unlike in all others sets of spin-dependent 
parton densities, is responsible for the opposite sign of 
dAa c /dy±dAa 5 /dy obtained with DNS (KKP). In the de- 
nominator of A Ac, gluon-gluon fusion does not drop out 
and can make a significant contribution depending on the 
amount of gluon polarization Ag. The set with the largest 
Ag, AAC, gives the largest cross section. The other three 
sets have relatively small gluon distributions and quark- 
antiquark annihilation is equally important. On the other 
hand, different sets of polarized parton densities have only 
very limited impact on the results for dAa c /dy±dAa c jdy 
and AAc for GSI-FAIR. This can be expected, since in pp 
collisions at small ^/S one predominantly probes the fairly 
well constrained valence quark distributions, and uncer- 
tainties in the polarized sea-quark and gluon densities do 
not matter much. 

Finally, we note that without taking into account any 
experimental cuts, Ac and AAc are anti-symmetric func- 
tions in the c.m.s. rapidity for pp collisions, as can be antic- 
ipated from Figs.[7J-[9](which, however, do include certain 



acceptance cuts). Similarly, Ac and AAc are symmetric 
in c.m.s. rapidity for pp experiments. This is a consequence 
of the anti-symmetric and symmetric initial states pp and 
pp, respectively. The charge asymmetry on the partonic 
level for the dominant qq subprocess implies that_ Q is 
preferentially emitted into the direction of q and Q into 
the direction of q. As was explained in [18j , in pp collisions 
one then finds an excess of centrally produced Q, while Q 
dominates at large absolute rapidities. This is also what we 
observe for J-PARC in Figs. [8j and EE The size of Ac and 
AAc for PP collisions at J-PARC is significantly smaller 
than for pp collisions at GSI-FAIR, simply because of the 
fact that for the relevant qq subprocess both partons are 
valence quarks in pp, greatly enhancing its relative con- 
tribution. When integrated over rapidity and without any 
kinematical restrictions, Ac and AAc vanish, and the to- 
tal yields of charm and anti-charm quarks are the same. 

3.4 Expectation for longitudinal spin asymmetries 

Experiments usually present their spin-dependent mea- 
surements in terms of spin-asymmetries rather than po- 
larized cross sections. The double-spin asymmetry All, 
defined in Eq. ([2J, has the advantage that many experi- 
mental uncertainties cancel in the ratio, in particular, it 
is not required to determine the absolute normalization 
of the helicity dependent cross sections. Also, theoretical 
uncertainties may cancel to some extent in All- However, 
before exploiting this, one has to make sure that pQCD is 
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Fig. 13. As in Fig.lTTIbut now as function of the c.m.s. rapidity 
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applicable in the relevant kinematical regime by compar- 
ing, for instance, the unpolarized cross section with data. 



Figure [TT] shows the integrated double-spin asymmetry 
A LL as function of p™ in at NLO accuracy for GSI-FAIR 
(left panel) and J-PARC (right panel), using the GRSV 
"standard" distributions and the same experimental cuts 
as for the underlying cross sections shown in Figs. [4] -[6] 
As usual, the scale uncertainty is shown for three differ- 
ent choices of m c , varying (i = fx r = (if in the range 
m™ < fi < 2m™ m . As can be seen, there is still a signif- 
icant scale ambiguity and, in case of J-PARC, also a de- 
pendence on m c . We refrain from showing LO results, but 
we note that because NLO corrections tend to be larger in 
the unpolarized discussed above, All is typically 

reduced by a factor of about two when NLO corrections 
are included. 

The sensitivity of All to different sets of polarized 
parton distributions is studied in Fig.[l2] As expected, the 
differences are small for GSI-FAIR which mainly probes 
the fairly well known valence distributions. Only at large 
p™ m , which corresponds to currently unexplored momen- 
tum fractions — > 1, some differences are noticeable. At 
J-PARC, expectations for All depend much more on the 
choice of spin-dependent parton densities, in line with the 
observations already made in the lower panel of Fig. [TO] 

Finally, in Fig. [13] we present the spin asymmetry as 
a function of the c.m.s. rapidity y and integrated over 
transverse momentum. The behavior of All for GSI-FAIR 
in the upper panel for positive y is driven by the cut on 
the charm momentum pi a b > 1 GeV. Both mass and scale 
uncertainties do not cancel and remain significant. 
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4 Conclusions 

We have performed a detailed study of the physics op- 
portunities with open charm production at low c.m.s. en- 
ergy pp and pp collisions at GSI-FAIR and J-PARC, re- 
spectively including unpolarized and polarized cross sec- 
tions and charge and spin asymmetries. All calculations 
are done at 0(a 3 a ) accuracy, and theoretical uncertainties 
due to the choice of scales fif >r and the charm mass m c 
are discussed in detail. In general, they turn out to be sig- 
nificant with the exception of the mass dependence of the 
charge asymmetries Ac and A Ac- 

Measurements of all these quantities would further our 
understanding of the perturbative QCD framework and 
the nucleon structure expressed in terms of unpolarized 
and polarized parton densities. The latter are probed at 
large momentum fractions x\_i, which are difficult to ac- 
cess at high-energy colliders. The charge asymmetry van- 
ishes at O{o? s ) and hence is a clear probe of non-trivial 
QCD dynamics beyond the LO. Detailed comparisons be- 
tween our theoretical expectations for unpolarized charm 
yields and future data will reveal to what extent perturba- 
tive methods are applicable in hadron-hadron collisions at 
VS w 10 15 GeV and will open up a window to pertur- 
bative resummations and/or the transition into the non- 
perturbative regime so far only little explored. 
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